Equilibrium and dynamics of a trapped superfiuid Fermi gas with unequal masses 
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Interacting Fermi gases with equal populations but unequal masses are investigated at zero tem- 
perature using local density approximation and the hydrodynamic theory of superfluids in the pres- 
ence of harmonic trapping. We derive the conditions of energetic stability of the superfiuid config- 
uration with respect to phase separation and the frequencies of the collective oscillations in terms 
of the mass ratio and the trapping frequencies of the two components. We discuss the behavior of 
the gas after the trapping potential of a single component is switched off and show that, near a 
Feshbach resonance, the released component can still remain trapped due to many-body interaction 
effects. Explicit predictions are presented for a mixture of 6 Li and 40 K with resonant interaction. 
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The superfiuid behavior of dilute interacting Fermi 
gases at very low temperature is now rather well under- 
stood both from the experimental and theoretical point 
of view [l| . In particular the attractive nature of the in- 
teraction between the two different spin components of 
the gas is known to play a crucial role along the whole 
BCS-BEC crossover. This includes the case of small and 
negative values of the scattering length, where the ordi- 
nary BCS regime of superfluidity holds, the BEC regime 
characterized by the formation of molecules in the Bose- 
Einstcin condensed state and the unitary regime where 
the scattering length takes a divergent value and the at- 
traction results in peculiar many-body effects. 

A more recent and intriguing direction is the search 
for superfluidity in mixtures of Fermi gases belonging 
to different species, and hence having different masses. 
Experimentally, the most promising candidates are ul- 
tracold mixtures of K and 6 Li, where the mass ratio 
is 6.7, near a heteronuclear s-wave Feshbach resonances 
. Equilibrium configurations of the uniform superfiuid 
phase where the atom densities of the two species are 
equal but the masses are different have been theoreti- 
cally investigated in Ref.0] within BCS mean field the- 
ory and, more recently, by quantum Monte Carlo meth- 
ods Other recent works have explored the interplay 
between different masses and different populations of the 
two components 0, Q . 

In this Letter we investigate the equilibrium and the 
dynamic properties of a trapped Fermi gas with unequal 
masses using local density approximation and developing 
the hydrodynamic theory of superfluids at zero tempera- 
ture. We assume that the external potentials for the two 
components, hereafter called f and J,, are harmonic and 
given by V£,(r) = m t (^i 2 + uj 2 y<y y 2 + u 2 za z 2 )/2 where 
a =t, I and mj and toj are the atomic masses. Since the 
two species have different magnetic and optical proper- 
ties, the trapping frequencies Ui a can be tuned separately. 



We consider a mixture of two different fermionic species 
with equal populations Nf = JVj. = N/2, correspond- 
ing to the most favorable condition for Cooper pairing. 
We assume that the gas is superfiuid at zero tempera- 
ture and we explore configurations where the densities 
of the two components are equal and move in phase: 
7if = = n/2, Vf = V| = v. We will not consider here 
exotic polarized phases, like the FFLO phase, whose rel- 
evance for trapped Fermi gases is still unclear at present. 

At equilibrium, where v = 0, the atomic density profile 
no(r) of the gas is given by the local density (also called 
Thomas- Fermi) approximation for the chemical potential 



Ho = n(n (r)) + V ho (r), 



(1) 



where /io is fixed by the normalization condition 
J no(r)dr = N and fi(n) = de/dn is the chemical po- 
tential of uniform matter, e(n) being the energy per unit 
volume of the homogeneous phase. In Eq.JT]) we have 
introduced the effective trapping potential 
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V ho (r) = |(V^ + V^y^l^ + ^ ya y 2 + Co^z'), (2) 
given by the average of the two potentials and 
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with m — (m-f + mjJ/2. In the superfiuid phase the 
densities of the two components are equal, even if the 
trapping potentials, in the absence of interactions, would 
give rise to different equilibrium profiles at zero tempera- 
ture, i.e. even if the oscillator lengths hj {m a u> a ) v l' 2 of the 
two components do not coincide. If the oscillator lengths 
are equal, the effective frequencies ^ simply reduce to 
the geometrical averages uji — y/UJj^LdjJ. 

Let us discuss the behavior of the equation of state 
along the BCS-BCS crossover and the corresponding 
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shape of the density profiles. At unitarity, where the 
scattering length a diverges, the equation of state takes 
the universal form 1] fi = (37r 2 ) 2 / 3 ft 2 (l + /?)n 2 / 3 /4m r , 
where m r — m^mj(m^ + m^) is the reduced mass and 
[3 is a dimensionless parameter depending on the mass 
ratio m-f/mj, and accounting for the interacting effects in 
this strongly interacting regime. By inserting /i = an 2 / 3 
inEq.©, with a = (37r 2 ) 2 / 3 (l+/3)/4m r , we find that the 
density distribution of the resonant gas takes the usual 
form 7i (r) = (po - Vho(r)) 3/2 /a 3/2 , the Thomas-Fermi 
radii Ri, where the density vanishes, being given by 



the new values 
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where Who = (w a; w !/ Wj) 1 '' 3 is the geometrical average 
of the three effective oscillator frequencies and = 
h/ Cj-h^m^mi) 1 / 2 . The value of (3 is known for the special 
case of equal masses m-f = mj., where (3 = —0.58 @, 0. 
Preliminary Monte Carlo calculations suggest that (3 de- 
pends very weakly on the mass ratio [9(. In the deep 
BCS limit, corresponding to a weakly attractive interac- 
tion (n|a| 3 -C 1), the equation of state and the Thomas- 
Fermi radii are given by the same expressions holding at 
unitarity by simply setting (3 = 0. 

In the limit of small and positive scattering length, 
corresponding to na 3 <C 1, the gas instead corresponds 
to a BEC of diatomic heteronuclear molecules of mass 
TOf + = 2m and molecular density n/2. In this 
regime the equation of state is fixed by the repulsive in- 
teraction between molecules and takes the usual bosonic 
form /j, m = g m n/2, where the coupling constant g m is 
related to the molecule-molecule scattering length a m 
by g m — 2irh 2 a m /m. The exact value of the molecu- 
lar scattering length has been calculated by Petrov et 
al. [l(| as a function of the atom scattering length a and 
the mass ratio m^/mp The bosonic chemical poten- 
tial is related to the fermionic one by fj, m = —Eb + 2/i, 
where Eb = —H 2 /2m r a 2 is the two-body binding en- 
ergy. From Eq.JT]), the density profile is then given by 
n(r) = (/iQ — Vh (r))2m/irh 2 a m corresponding to the 
Thomas Fermi radii 
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where il>h = {Cj x Cby(jj z ) 1 / 3 and a^ c = H/ yJ2muj\ lo . No- 
tice that a^ BC differs from the oscillator length aho de- 
fined above for the resonant case. 

We now take advantage of the fact that the trapping 
potentials of the two different species can be tuned sepa- 
rately to suggest an experiment pointing out in a direct 
way the attractive role of the interactions in the pres- 
ence of a Feshbach resonance. After generating the equi- 
librium configuration discussed above we switch off the 
confining potential of a single species, say Vj* = 0, corre- 
sponding to a change of the trapping frequencies ([3]) into 



TTlf + 771J, 



(6) 



The potential can be switched off either adiabatically, 
bringing the system into a new equilibrium configuration, 
or suddenly, giving rise to the excitation of collective os- 
cillations (see discussion in the second part of this work) . 

In the absence of interactions, the |-atoms would fly 
away leaving an ideal gas of |-fermions trapped in the 
harmonic potential VJ . This will be also the case in 
the deep BCS superfluid regime where interactions are 
too weak to keep the j-fermions confined. In the other 
(more robust) superfluid regimes, however, the released 
atoms do not necessarily escape to infinity but can remain 
trapped due to the attractive interaction with the other 
species. This statement is obvious in the BEC regime, 
where each |-fermion forms a bound molecule with a 
corresponding |-particle. At unitarity, however, the two- 
body binding energy Eb = —h 2 /2m r a 2 vanishes meaning 
that no molecules can exist in vacuum. In this case the 
trapping of the J, component is a pure many-body effect 
reflecting the attractive nature of the interatomic force. 

It is not difficult to derive explicit conditions for the 
energetic stability of the new superfluid configuration. 
The stability is ensured if the energy Eg of the config- 
uration where the two components remain trapped and 
fully overlapped is smaller than the energy En of the 
normal state where the gas is phase separated and only 
the T-atoms are trapped. At unitarity the energy of the 
trapped superfluid state is given by 
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and, in the absence of trapping for the |-atoms, one 
has cD/j = yj m-f / (m-[ + m^^x^LOy^uj^) 1 ^ . Since the 
energy of a trapped gas of non-interacting ] fermions 
is given by En = h{uj x ^ujy^u} z -\ ) 1 / 3 (3A r ) 4 / 3 /8, we find 
that the superfluid gas, where both components remain 
trapped, is energetically stable if the condition 
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is satisfied. Taking into account that (3 barely depends 
on the mass ratio and hence remains close to the 
equal mass value (3 = —0.58, the above condition is 
always satisfied if m| < mj., showing that the super- 
fluid remains energetically stable if we release the po- 
tential of the heavy species. Conversely, if the condi- 
tion ||5J) is violated, the superfluid configuration corre- 
sponds to a metastable state which is energetically un- 
stable toward phase separation. It is also interesting to 
compare the Thomas Fermi radii Ris and Rin of the 
trapped cloud in the (new) superfluid and in the sepa- 
rated normal phase, respectively. A simple calculation 
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yields R iS / R iN = (l + /?) 1/4 ((m T + m i )/m i ) 1 / 4 , showing 
that the superfluid phase corresponds to the configura- 
tion with smaller radii if and only if Eq.© is satisfied. 

Let us now discuss the macroscopic dynamic behav- 
ior of the superfluid. This is is obtained by deriving 
the hydrodynamic equations of motion in terms of the 
atom density n(r,t) and the velocity field v(r,i). The 
Lagrangian L of the system, in the local density approx- 
imation (LDA), is given by 
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where (f> is the phase of the order parameter 
(* T (r,i)* i (r,i)) = K^^r)*;^))^^*), * ff (r,t) being 
the fermionic field operators. The superfluid velocity v 
and the phase <j) entering the Lagrangian @ are related 
by the most important condition 
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(10) 



where m-f + toj, is the mass of the pair, which can be 
derived microscopically by noticing the state moving with 
velocity v is obtained from the steady state by applying 
the gauge transformation y^fr) — > \& a (r)e m,yV ' r ^ h to the 
two Fermi field operators 

Taking Eq. (fTU|) into account, the equations of motion 
of Lagrangian ([9]) yield the hydrodynamic equations 
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+ V(nv) = 0, 



d~v f 1 



(11) 
0, . (12) 



Equations (fTTj) and (|12p apply to weakly as well as 
to strongly interacting superfluids and permit to calcu- 
late the macroscopic dynamics of the system (expansion 
and collective oscillations) once the equation of state is 
known. 

In uniform matter (Vh = 0) the linearized solu- 
tions are phonons with sound velocity fixed by c 2 = 
(n/m)dfj,/dn. At unitarity, the sound velocity is given 
by c = ftfc Fv /l/(3m T TOjVl + /?, where k F = (S^n) 1 ^. 
In the BCS limit, the same formula holds with (3 = 

For trapped configurations, the frequencies of the col- 
lective oscillations are obtained 13] by linearizing the 
hydrodynamic equations around the equilibrium distri- 
bution no(r) and depend on the applied external poten- 
tials only through the effective oscillator frequencies ((3]). 
They also depend on the equation of state and in partic- 
ular on its density dependence [14|]. A special case is the 
center of mass oscillations of the cloud whose frequency, 
in the superfluid phase, is independent of the equation of 
state and given, for each direction, by u>i. We emphasize 
that in general tDj differ from either u>^ and u>ii, point- 
ing out that in the superfluid phase the two components 
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FIG. 1: Time evolution of the size of the cloud along the 
radial direction after switching off suddenly the trapping po- 
tential for 40 K fermions. Here r = £)±(Q)t. The dashed line 
corresponds to the new equilibrium value (see text). 



always oscillate in phase as a consequence of the pairing 
mechanism. 

Equations (fTT|) and (112|) remain valid even if the effec- 
tive trapping frequencies depend explicitly on time, i.e. if 
LUi = LUi(t). Remarkably, if the chemical potential has the 
power law dependence /i oc n 7 on the density and the con- 
fining potential is harmonic, the hydrodynamic equations 
admit a class of exact solutions given by v(i, r) = at-r and 
n(t, r) = n (x/b x ,y/b yi z/b z )/ (b x b y b z ), where o,(i), 
are time-dependent parameters. From the scaling form of 
the density, we see that the parameters 0i(t) are related 
to the Thomas Fermi radii Ri(t) of the evolving cloud 
according to Ri(t) — Ri(0)bi(t). Inserting the scaling 
ansatz in Eq. (|ll|) yields the relationships Ui — bijbi and 
one then obtains a set of coupled ordinary differential 
equations [lit 



u>{(t)bi 



^ 2 (0) 



1 



bi (b x b y b z p 



(13) 



Equations (|13[) apply both to resonance, where 7 = 2/3, 
and to the BEC regime where 7 = 1. 

If one suddenly switches off the trapping potential for 
the ^-component, corresponding to suddenly setting the 
effective frequencies to the new values © , the superfluid 
gas will start oscillating around the new equilibrium con- 
figuration. The dynamics of the Thomas Fermi radii of 
the cloud can be calculated from Eq. (fTB"]) by substituting 
u)i(t > 0) = cl>mcw and employing the initial conditions 
&i(0) = 1 and bi(0) = 0. 

For simplicity, let us assume that the two trapping po- 
tentials are axis-symmetric (io X o = u ya = u)± a ). We fur- 
ther assume that the two components have initially equal 
oscillator lengths (and hence, in the degenerate limit, 
equal density profiles also in the absence of interaction) , 
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FIG. 2: Time evolution of the size of the cloud along the 
radial direction after switching off suddenly the trapping po- 
tential for 6 Li fermions. Here r = u>±(0)t. The dashed line 
corresponds to the new equilibrium value (see text). 



configuration is no longer symmetric. We emphasize that 
in this second case the superfluid configuration is ener- 
getically unstable. A major question in this case is to 
understand the decay mechanisms and the role played 
by the sudden excitation of the collective modes. 

In conclusions, we have derived the equilibrium condi- 
tions and hydrodynamic equations of a superfluid Fermi 
gas with unequal masses and investigated the behavior of 
a gas at unitarity after the release of the trapping poten- 
tial of a single component. We have shown that, under 
appropriate conditions, the superfluid phase is energeti- 
cally stable against phase separation of the two compo- 
nents. As a result, the released fermions remain confined 
in the trap due to the pairing with the other compo- 
nent, pointing out in a direct and remarkable way the 
attractive nature of the interatomic forces near a Fesh- 
bach resonance. 

This work is supported by the Marie Curie program 
under contract EDUG-038970 and by the Ministero 
dell'Istruzione, dell'Universita e della Ricerca (M.I.U.R.). 



b±, Eqs (1131 take the form 



corresponding to o>^(0) = oj±-\uj±i and Co z {0) = uj z -\uj z i. 
By introducing the dimensionless time r = uj±(0)t and 
setting b x = b 
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where rj = m,[/(m^ + mjj and A = u> z {0)/Q±(0) is the 
aspect ratio of the trapping potential (|2|) . 

We have solved Equations (fl"4|) and (fl~5f for a mix- 
ture of 40 K and 6 Li fermions in an elongated trap with 
A = 0.2. In Fig[T] we plot the calculated time evolu- 
tion of the radius of the cloud in the radial direction, 
after we suddenly switch off the trapping potential of 
40 K fermions, corresponding to rj = 0.869. The static so- 
lution b± = b z — 1/n 1 / 4 of these equations and the cor- 
responding Thomas- Fermi radii Ri ncw — i?i(0)/?7 1 / 4 = 
1.036i?i(0) characterize the new equilibrium conditions 
after releasing the J. trapping potential. Since rj is close 
to 1, the initial configuration is close to equilibrium 
and the resulting oscillation is a linear superposition of 
the two breathing modes, in the radial and longitudi- 
nal directions. For elongated cloud, corresponding to 
A <C 1, these are given by w ra d = ^/10/37y 1 / 2 (Dj_(0) and 

^axial = yi27V/ 2 A^(0) fl. 

In Fig[2] we instead plot the time dependence of the 
transverse radii of the cloud, after we suddenly switch off 
the trapping potential of the 6 Li fermions, corresponding 
to rj = 0.131 and Ri ne w = 1.66i?i(0). In this case, the 
initial configuration is far from equilibrium and non lin- 
ear effects play an important role. In particular we see 
that the breathing of the cloud around the equilibrium 
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